Electro-hydrodynamics of binary electrolytes driven by modulated surface potentials 
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We study the electro-hydrodynamics of the Debye screening layer that arises in an aqueous binary 
solution near a planar insulating wall when applying a spatially modulated AC-voltage. Combining 
this with first order perturbation theory we establish the governing equations for the full non- 
equilibrium problem and obtain analytic solutions in the bulk for the pressure and velocity fields of 
the electrolyte and for the electric potential. We find good agreement between the numerics of the 
full problem and the analytics of the linear theory. Our work provides the theoretical foundations 
of circuit models discussed in the literature. The non-equilibrium approach also reveals unexpected 
high-frequency dynamics not predicted by circuit models. 

PACS numbers: 47.65.+a, 47.32.-y, 47.70.-n, 85.90.+h 
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I. INTRODUCTION 

Recently, there has been quite some interest in electro- 
hydrodynamics in microfluidic systems. AC-driven, mod- 
ulated surface potentials have been used for pumping, 
fluid circulation, and mixing [1-13]. For an overview of 
AC electro-osmosis we refer to Refs. 14-19 and references 
therein. 

Although AC electro-osmosis is typically analyzed with 
the help of homogeneous circuit elements open questions 
remain about the applicability of such approaches [19]. 
We revisit the problem studied by Ajdari [2] where an 
electrolyte is perturbed by an AC-driven spatially mod- 
ulated surface potential, but include explicitly an insu- 
lating layer between the electrode providing the driv- 
ing potential and the electrolyte. We primarily think 
of this insulator as an oxide grown intentionally for de- 
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FIG. 1: A sketch of the system under study. The binary 
electrolyte is situated in the half space x > 0. Below it, 
for — d < x < 0, is a planar wall consisting of an insulating 
dielectric slab of thickness d and below that, for x < —d, 
is a semi-infinite conductor. The top surface, x = — d, of 
the conductor is biased by a periodically modulated potential 
V ext (y, t) of period 2n/q (dotted line), which gives rise to the 
formation of a Debye screening layer of thickness \ D in the 
electrolyte (dashed line). 



vice purposes, but it could also represent the molecular 
Stern layer in case of non-oxidized electrodes. We de- 
velop a full non-equilibrium description of the electro- 
hydrodynamics of this system thus extending previous 
modeling of the surface and the Debye layer as simple 
capacitors. This allows us to study the full dynamics of 
ion concentrations, electrical potentials, velocity fields, 
pressure gradient-fields, and electrical currents as well as 
the justifications for a description based on homogeneous 
circuit elements. 

In the following we consider a binary electrolyte, i.e., 
an aqueous solution of a salt containing a positive and 
a negative type of ions with charges +Ze and — Ze, re- 
spectively, where Z is the valence and e the elementary 
charge. In terms of Cartesian coordinates xyz the elec- 
trolyte is confined to the semi-infinite space x > by an 
impenetrable, homogeneous and planar insulating layer 
with dielectric constant e s placed at —d < x < 0, see 
Fig. 1. The insulating layer is bounded by a conductor 
at x < — d which has been biased at the surface x = — d 
by a spatially modulated, external AC potential V^, xt (y, t) 



V tsst (y,t) = V cos(qyy 



(1) 



where Vq is the amplitude, q the wavenumber of the spa- 
tial modulation, and u> the driving angular frequency. 

There is complete translation invariance along the z 
axis, so the z coordinate drops out of our analysis, and 
all positions r = xe x + ye y are therefore just referring to 
the xy plane. 

The manuscript is organized as follows: in Sec. II we 
present the non-equilibrium description and in Sees. Ill 
and IV we analytically study linearized equations of the 
static and dynamic regimes, respectively. In Sec. V we 
study numerical solutions of the fully coupled non- linear 
electro-hydrodynamic problem. Finally, in Sec. VI we 
compare these solutions with the analytical solutions of 
the linearized equations, and furthermore contrast our 
results with the literature, before we in Sec. VII conclude. 



2 



II. NON-EQUILIBRIUM DESCRIPTION 

The basic non-equilibrium formalism for continuum 
electro- hydrodynamics is well-known (see, e.g., Ref. 19), 
but as mentioned in the introduction we explicitly in- 
clude an insulating layer in the description. We do not 
include an intrinsic zeta-potential, i.e., no un-passivated 
surface charges on the insulator-electrolyte interface. We 
note that experimentally any intrinsic homogeneous zeta- 
potential may be compensated by a corresponding DC- 
shift added to the applied AC potential. Due to hetero- 
geneous surfaces, it may be anticipated that although the 
average zeta-potential is nulled out, there might be fluc- 
tuations left. These will be the topic for future work. For 
zero intrinsic zeta-potential we solve the full non-linear 
equations numerically, but to obtain analytical result we 
also study the linearized equation with special emphasis 
on the capacitance due to the insulating layer. In Sec. V 
we show that for experimental relevant parameters, the 
linear theory is surprisingly good. Typical values of var- 
ious central parameters are listed in Table I. 

A. The insulating layer, — d < x < 

The insulating layer contains neither free space charge 
nor free currents so the electrical potential </>(r, t) is gov- 
erned by the Laplace equation, 

V 2 0(r, t) = 0, for - d < x < 0. (2) 

B. The electrolyte, x > 

In the liquid electrolyte we consider the ionic densities 
n ± (r,t), the potential <f>(r, t), the ionic current densities 
(the ionic flux densities) i ± (r, t), the velocity field v(r, t) 
of the electrolyte, and the pressure p(r, t). In the follow- 
ing we suppress (r, t) unless needed for clarity. 

The number densities of the ions couple to the poten- 
tial via Poisson's equation, 

V 2 0=-^(n+ -n"). (3a) 

The ionic current densities are coupled to the ionic 
densities by a continuity equation, which in the absence 
of any chemical reactions in the system is 

a t n ± =-V-i ± . (3b) 

The presence of convection or of gradients in the den- 
sities and the electric potential <f> will generate ionic 
current densities i*. The Nernst-Planck equation gives 
these currents 

i^-BV^+n^Tfi^Vci, (3c) 

where, for simplicity, we have assumed that the two types 
of ions have the same diffusivity D and the same mobility 



TABLE I: Typical values of central parameters. 
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Debye length 
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Resonance frequency 
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Debye frequency 
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rad/s 


Critical frequency 
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Convective voltage 


T/ 
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Ionic density 
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molL" 1 


Viscosity 
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KT 3 


Pas 


Mass density 




P 


10 3 


kgm" 3 


Ionic diffusivity 




D 


KT 9 




Capacitance ratio 
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\x. We remind the reader that both the diffusivity D and 
the electric conductivity a are linked to the mobility [i via 
the Einstein relation D = (k B T/Ze)Li and a ± = Zen^Li. 

Finally, the velocity field and pressure of the liquid 
are coupled to the potential and ionic densities by the 
Navier-Stokes equation 

p[d t v + (v • V)v] = -Vp + r/V 2 v - Ze[n + - n~] V<j), 

(3d) 

where p is the mass density, r\ is the viscosity of the 
liquid, and p is the pressure. Furthermore, treating the 
electrolyte as an incompressible fluid we have 

V • v = 0. (3e) 

The coupled field-equations, Eqs. (3a) to (3e), fully gov- 
ern the physical fields n ± , 0, i ± , v, and p. 



C. Boundary conditions 

Assuming a vanishing zeta-potential the boundary con- 
dition for the electric potential is 

<KT,t)\ x= _ d = V ast {y,t), (4a) 
^M)L =oo = 0. (4b) 

At the interface between the electrolyte and the insu- 
lating region the normal component of the ionic current 
density vanishes, 

= a x n± (r, t) | x=0 ± (r, i)0 a 0(r, t) \ x=Q . (5) 

Here, we have utilized Eq. (3c) and the absence of con- 
vection at the interface due to the no-slip boundary con- 
dition, 

v M)Lo = °- ( fi ) 
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For the ionic densities we have 



n (r, t)\ = n 

y ' '\x=oo ° 



(7) 



where is the homogeneous density of either of the two 
types of ions in the absence of an external perturbation, 
i.e., when Vo = 0. For the pressure, we assume that 
we have no externally applied pressure gradients so that 
p is the internal pressure caused by fluid flow and the 
electrical forces on the ions. 



III. STATIC REGIME, w = 

In the static regime we have equilibrium and neither 
current nor fluid flow, i.e., i =0 and v = 0. The pres- 
sure gradient balances the electrical forces on the charges. 
The governing equations for <j> and of course reduce 
to those of electro-statics. 

In the insulating layer it follows from Eqs. (2) and (4a) 
that 



0(r) 



B ie - qx +B 2 e qx 



cos(qy), for— d<x<0, (8) 



where B l and B 2 are integration constants. 

In the electrolyte 0(r) is governed by the nonlinear 
Poisson-Boltzmann equation [19] 



V 2 0(r) 



k T 
B sinh 



Ze\ 2 D 



Ze 



for x > 0, (9) 



expression relating the total potential drop across the 
system and the accumulated charge in the electrolyte, 

roc 

0(00, y) - <K-d,y) = C-g 1 / dxZe[n+(v) - n"(r)]. 

Jo 

(13a) 

The coefficient, 

^ = [1 + ^^^ 



+ y/l + (q\ D ) 2 COsh{qd) C D \ (13b) 

is identified as the inverse of an effective series capaci- 
tance. The constant C s is the intrinsic surface capaci- 
tance and Cd the capacitance of the Debye layer given 
by 



r = — 
a 

C D .f. 
Ad 



(13c) 
(13d) 



In Ref. 2 the potential in the bulk of the electrolyte 
(x 3> Ac) is governed by the Laplace equation, which 
is coupled to the external potential V ext by an effective 
capacitance C given by 



Co^iC-' + Cn 1 )- 1 . 



(13e) 



introducing the Debye screening length 



It follows from Eq. (13b) that this approach for u = is 
valid up to second order in the parameters q\o and qd. 



IV. LINEARIZED DYNAMIC REGIME, w > 



(10) 



( 2Z 2 e 2 nJ 

For q going to zero we have a constant surface potential 
(r)L ^</>o (Ha) 



lim ( 

q^0 



and the solution to Eq. (9) is given by the well-known 
Gouy-Chapman solution [19] 



lim c/)(r) 

q^0 



Ak^T 
Ze 



arctanh 



tanh 



4fc B T 



-x/X D 



(lib) 

For g^Owe are not aware of any analytical solutions, 
but as we shall show, analytical results can be obtained 
in the Debye-Hiickel approximation Ze<f> <C k B T, where 
Eq. (9) becomes linear, 



(r) = A, 



for x > 0. 



(12) 



Here, the corrections are to third order in Ze<j>/k B T be- 
cause sinh(a;) = x+C(x 3 ) . The space charge follows from 
Poisson's equation, Eq. (3a). From a straightforward so- 
lution for (j) and Ze(n + — n~) we arrive at the following 



We now solve Eqs. (3) in the dynamic regime, oj > 0. 
First the ionic current densities are eliminated by insert- 
ing Eq. (3c) into Eq. (3b). Using the incomprcssibility of 
the fluid, Eq. (3e), we get the continuity equation 

= D V 2 ^ - (Vn ± ) ■ v ± /iV • (n ± V<j>). (14) 



A. Debye— Hiickel approximation 

To advance further by analytical methods, we now lin- 
earize the continuity equation, Eq. (14), in the density as 



follows. We assume n (r, t) _^ 



and write 



n ± ( r , t) = + < 5n ± (r, t), lim (5n ± (r, t) = 0. (15) 

Since we assume a zero intrinsic zeta-potential it is a 
non-zero Vq that spawns Sn ^ 0, and when the applied 
voltage Vo is much smaller than the thermal voltage Vt, 
defined by Vt = (1 + CnlC s )k B T jZe (as we shall see in 
the next subsection), we have \&rv si \ <C n x . In this limit 
the Debye-Hiickel approximation is valid, and V0 is 
substituted by n ro V(/> in Eq. (14). We subsequently use 
Eq. (3a) to replace V 2 with —Zev/e where 



v = n 



5n + — 5n 



(16) 
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Finally, we form the difference of the "±" -versions of 
Eq. (14) and obtain the partial differential equation 



d f u = 



DV 2 -%-v.V 



V. 



(17) 



Together with the boundary condition for (p at x = — d 
as well as the continuity of 4> and ed x (f> at x — we 
may determine the constants straightforwardly keeping 
in mind that k depends on q and u). For the coefficients 
in the electrolyte we get 



B. Diffusive regime 



Ci = -= — C eS (u), 

Ze q 



(23a) 



Our study of the static regime reconfirm the well- 
known result that the net charge density is non-zero only 
in the Debye layer, x < 3Ad. In this region convection 
will be suppressed due to the no-slip boundary condi- 
tion. Thus, convection can be neglected, diffusion will 
dominate (corresponding to a low Peclet number), and 
the electrodynamics can be solved independently of the 
hydrodynamics. On the other hand, the hydrodynam- 
ics of course still depends on the electrodynamics via 
the body- force. Since the density difference v changes 
over the length scales Xd and q^ 1 for the x and y di- 
rections, respectively, the condition for the decoupling is 
l^xl/A-D + \v v \q <C Dq 2 for < x < 3X D . In this limit 
Eq. (17) has a general cos(qy)e luJt modulated decaying 
solution of the form 



and 



v = de KX cos(qy)e i 



x>0, 



(18a) 



where the decay parameter k depends on the ratio be- 
tween the frequency uu and the Debye frequency u D , 



K = T^" Jl + (q\ D ) 2 +i—, 

D 



(18b) 
(18c) 



For the potential we seek a solution of a form similar to 
Eq. (18a), <f> oc cos{qy)e lu>t , and substituting this together 
with Eq. (18a) into Eq. (3a) yields 

{d 2 x - q 2 )<P = -^C ie - ra cos(^)e-*. (19) 

Demanding </>(r, t)\ = the solution is 

Ze/e 



q 2 — k 2 



C 2 e- qx 



cos(qy)e 1 ' 



x>0. 



(20) 

In the insulating layer we have the following cos{qy)e lu ' t 
modulated general solution to Eq. (2), 

4> = \c 3 e- qx + C 4 eH cos(qy)e i0Jt , -d < x < 0. (21) 

In order to determine C\, C2, C3, and C4 we first con- 
sider the boundary condition for the current. Applying 
the Debye-Hiickel approximation to the second term in 
Eq. (5) and forming the difference of the " ±" solutions 
we arrive at 



= d x 



i/(r,t) + 



C. 



D 



Ze\ D 



ct>{v,t) 



x=0 



(22) 



C 2 = -1 Ci, 

u D q 



(23b) 



while for the insulator the coefficients have a similar, but 
less compact form. Above, 



! , 2 sinh(grf) x 

6 off ( W ) = ( kA -D) ~ d C s 



(23c) 



q\ D {q\ D + k\ d ) + i ^f- 

+ k\d r— r-r ; , - cosh(qd) C D , 

q\ D {q\ D + n\ D ) 

which satisfies the definition in Eq. (13a) and reduces to 
Eq. (13b) in the DC limit. From the general solution for 
the potential in the electrolyte, Eq. (20), we may now in 
more detail examine the constraints on Vo for the Debye- 
Hiickel approximation to be valid. Straightforward cal- 
culations show that max{</>} <C k B T/Ze corresponds to 
V < V T = (1 + C D /C s )k B T/Ze for qd < 1, q\ D < 1, 
and low frequencies. 



C. Long-period and low- frequency modulation 

Next, we consider the regime where the spatial period 
of the modulation is much longer than all other length 
scales, i.e., q\o <C 1 and qd -c 1. We also assume that 
iv <C uo D so that k ~ 1/Xd- In this limit we get 



Ze\]j u* + iu: 



c/\ D 



cos(qy)e wt + 0([q\ D } 2 ), 



and 



• = V - 



iu> 



-e~' lx cos(qy)e^ t + 0(q\ D ) 



UJ* + VjJ 

where we have used the notation of Ajdari [2] 



resonance frequency: lu* = q\ D (l + 5) lu d , 
conductivity: a, 

capacitance ratio: 



(24) 
(25) 

(26a) 



6= c7- 



euj D , (26b) 
(26c) 



These results are equivalent to those in Rcf. 2 if we 
similarly to Eq. (13a) introduce the Debye layer surface 
charge a D (y) = Ze J °° dx u(x, y). 
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D. Body- force 



E. Linearized flow in quasi-steady state 



Until this point we have used the exponential nota- 
tion for the temporal dependence. However, since the 
body-force is essentially non-linear in the electrical po- 
tential/density [see last term in Eq. (3d)] we have to take 
the real part to get the body-force, i.e., F = —Zev'Vcf) = 
— ZeRe{^}Re{V0} so that we get 



rjv 1 cos(2wi + tp) 



-x/\d 



(27a) 



x [2cos 2 (qy)e x + sin(2qy)e y ] +0([q\ D ] 
where following Rcf. 2 we have introduced 



Vi EE 



4r?(l + 5) 

and the frequency dependent phase shift 



(27b) 



tp- 



arctan 



UJ 

2u? 



-) 

2uj J 



In the derivation of Eq. (27a) we have used that 



Re 



Re 



— cos(2ujt + tp) 
2lu*(^ + ' 



(28) 



At low frequencies, F oc w, it becomes maximal at the 
resonance frequency u>*, and then it falls off again at 
higher frequencies. We note that lim^^o F = 0([q\D] 2 ) , 
but this small force will just be balanced by a pressure 
gradient so that lim w ^ v = and lim w ^ = 0. 



In order to solve the Navier-Stokes equation, Eq. (3d), 
we note that for a body-force of small magnitude and 
with slow temporal variation the fluid response is lin- 
ear and the flow will approximately be at steady state at 
each moment in time. We begin by comparing the inertial 
terms on the left-hand side (LHS) with the viscous term 
(second term) on the right-hand side (RHS). The body 
force has a characteristic frequency u> and two character- 
istic length scales Ad and q~ 1 for the x and y-directions, 
respectively. Since dt essentially gives a factor of u>, and 
V essentially gives e x +qe y , we can show that the vis- 
cous term dominates over the LHS when ui <C uo c where 



uj c ee - mm{q , A D }• 
P 



(29) 



(27c) p or q\ D this means that uj c 
small Reynolds numbers, we get 



^ q 2 . In this way, for 



= -Vp + r?V 2 v + F, w<w c 



(30) 



which is the resulting quasi-steady flow problem which 
is linear in the velocity field. Slip-velocity approaches 
usually rely on this equation, see Rcf. 2 and references 
therein. However, Eqs. (30) and (27a) can actually be 
solved exactly with a solution given by 



v(r,i) = vi 



Vl 



cos(2wt + tp) 

UJ I UJ* 
UJ* UJ 

cos(2wt + tp) 

UJ I UJ^_ 

UJ* UJ 



-2qx 



2qX D G 1 (X D 1 ) cos(2qy)e x + Qi(2q) sin(2qy)e y 



(31a) 



2qxe 2qx cos(2qy)e x + {(1 - 2qx)e 



2qx 



l/\o ' 



}sm(2qy)e y +0(q\ D ) (31b) 



P (r, t) = -^ V rT t+ J ] ^ cob(2 W ) + 0,(0)] , 



UJ I UJ^_ 

UJ* UJ 



(l + 2qX D ) 2 



.^ COS^t + tp) ( e -2 qXcos{2qy)+ ^_^ [l+cos{2qy)] \ +0{qXDl 



UJ _|_ uf_ 

UJ* UJ 



4qX D 



(32a) 
(32b) 



as may be verified by direct insertion. Above, 
1 - 2q\ D 



Gi(k) = 



and 



4(?A D 



(33b) 



[i - (2q\D) 2 r 

x A e -(\~ 1 -2q)x q 2 \ )kx\ (33a) have been introduced. As seen the flow decays exponen- 
V ' / tially over a length scale of l/q when uj <C lo c . When 
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FIG. 2: The potential (j>, pressure p, and velocity field v. Panel (a) shows a gray scale plot of the amplitude of the potential <j> as 
a function of qx and qy, Eq. (25), and panel (b) the pressure p, Eq. (39). Notice the period doubling in the pressure compared 
to the electric potential. Panel (c) shows a snapshot of the harmonically oscillating velocity field v in the bulk, Eq. (38), and 
panel (d) likewise in the Debye layer, Eq. (35). The flow pattern contains rolls, which are indicated by contours of constant 
velocity (dashed lines). 



the frequency becomes comparable to or larger than ui c 
we have competing length scales since the <9tv term in- 
troduces an additional length scale, (rj/puj) 1 / 2 , which as 
mentioned becomes (uJc/uj) 1 ^ 2 x l/q for q\o -C 1. So 
in the above expressions for v and p we expect that the 
spatial cut-off length changes from (2q)~ 1 to A with 



A(w) — minjl, \fu c /u j. 



(34) 



Even for uj <C uj* the condition uj <C (jJ c is not necessarily 
fulfilled. In fact, for the numbers in Tabl e I we have 
uj c < uj* < u>d so at resonance 2qh{uj*) = ^Juj c /uj* <C 1. 

F. Flow and separation of length scales 

As mentioned above the flow is typically analyzed by 
slip-velocity approaches and here we show how such an 
approach gives an asymptotic solution in full agreement 
with the exact solution. We study the flow over a Afl- 
scale at the boundary first and then a q _1 -scale. For this 
boundary layer approach, we assume that for x < 3Ad, 
we have v x rs 0. Solving for the pressure and substituting 
into the y-component of Eq. (30) we get 

v v =v s (y,t) (l-e- x / x °)+0(q\ D ), x<3\ D (35) 
with the prefactor 

cos(2 UJ t + ip) . . . 
v s (y,t) ee Vl ^ ^ sin(2gy). (36) 

In the limit x > 3Ad and qXjj <C 1 the velocity v s can 
be interpreted as a slip-velocity at the wall acting as a 
conveyor belt for the bulk fluid, see Fig. 2(d). 



For x > 3Ac we have that F is exponentially sup- 
pressed and we solve Eq. (30) together with Eq. (3e) and 
the boundary condition 



V M)L =0 = v s(y,t) e y - 
To lowest order in qXu this gives 
cos(2 ujt + ip) ^_ 2qx 



(37) 



(38) 



^ — 2qx cos(2qy)e x + (1 — 2qx) sin(2qy)e y 



and 



p ~ —Aqrjv 1 



cos(2wt + (p) 



-2qx 



cos{2qy). 



(39) 



If we now substitute into Eq. (3d) we get RHS — 
LHS oc e ~ x ' XD + 0(uj/uj d ) + 0{[q\ D ] 2 ) which shows 
that Eqs. (38) and (39) are indeed excellent approxima- 
tions to the full solution of the non-linear time-dependent 
Navier-Stokes equation, Eq. (3d), for x > 3 Ad- For the 
incompressibility constraint, Eq. (3e), our solution gives 
V • v = 0([qX D ] 2 ). In Fig. 2(c) we show a plot of the 
velocity-field, Eq. (38), along with the contours for con- 
stant velocity. 

We note that in the limit x > 3A D the exact solutions, 
Eqs. (31b) and (32b) reduce to Eqs. (38) and (39) for the 
bulk. 



V. FULL NUMERICAL SOLUTION 

In this section we present results from numerical finite 
element simulations (Femlab) of the five coupled equa- 
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FIG. 3: The real and the imaginary parts of the normalized 
effective capacitance C e g(Vo,u))/Co versus the normalized fre- 
quency lu/uj* for four different values of the voltage amplitude 
Vo/V T . The dashed lines show the analytical results of the 
linear theory, Eq. (23c). The inset shows the relative de- 
viation (below resonance) of the full numerical solution for 
C e ff (Vo) /C from the Debye-Hiickel result as a function of 
V /Vt. 



tions, Eqs. (3a) to (3c), with the boundary conditions in 
Eqs. (4a) to (7). For simplicity we assume a low Reynolds 
number so that we can neglect the inertial term, (v • V)v, 
in Eq. (3d) . This provides full access to the temporal and 
spatial evolution of the physical quantities n ± 1 <p, i^, v, 
and p. The spatially periodic problem (in the y-direction) 
would typically be handled by applying periodic bound- 
ary conditions to the unit-cell (e.g., < y < 2n/q). How- 
ever, due to the symmetry of V ex t the computational do- 
main can be reduced to0<y<7r/g with homogeneous 
Dirichlet or Neumann boundary conditions. For the x di- 
rection our domain is cut off at a distance x = 6ir/q from 
the interface using Dirichlet boundary conditions for the 
fields. Near the interface to the insulator, < x < 3Arj, 
we employ a structured grid to resolve the Debye-layer. 
For the temporal evolution we employ the Femlab time- 
stepper directly starting from initial solutions at t = 
which are zero everywhere. The duration of the transient 
depends on the inertia in the system, but typically the 
temporal harmonic state is fully evolved after a time of 
the order 10 — 100 x uj^ 1 . For our simulations we have 
used the typical values in Table I except for the modula- 
tion where we have used q^ 1 = 10~ 6 m and consequently 
uj* = ui c = 10 6 rad/s. For the external potential V cx t we 
have used the imaginary part of Eq. (1) [which is zero at 
t = 0] with an amplitude < V < 10 V and a driving 
frequency 10 4 rad/s < uj < 10 8 rad/s. 

In order to directly compare our numerical results to 
the linearized theory we normalize frequencies with the 
resonance frequency uj* , velocities by v\, voltages by the 
thermal voltage Vt, and capacitances by Cq. 



Figure 3 shows numerical results for the effective ca- 
pacitance C e ff, see definition in Eq. (13a), as a function 
of the frequency for varying amplitudes Vq of the exter- 
nal voltage. The dashed line shows the corresponding 
analytical result from Eq. (23c). As seen there is a good 
agreement between numerical results and our analytical 
predictions even for Vq > Vt where the Debye-Hiickel 
approximation is typically expected to work poorly. We 
furthermore note that at low frequencies Re{C c ff} ap- 
proaches Co in full agreement with the analytics and the 
log- log plot also reveals two distinct regimes for uj < uj* 
and uj > uj* . In fact, dissipation is maximal exactly at 
the resonance frequency uj* predicted by the linear the- 
ory. The inset shows the relative error of the Debye- 
Hiickel approximation which at large voltages saturates 
at a value of the order (1 + S)^ 1 , here equal to 0.09. 

The linear theory predicts a harmonic velocity field 
with a vanishing time average and our numerical 
simulations confirm this low-frequency dynamics, see 
max r t {v(r, i)} in Fig. 4(a). The corresponding solid 
line shows exact results within the Debye-Hiickel ap- 
proximation [20] and the dashed line shows (-^ + 
as suggested by Eq. (38). As expected the induced 
harmonic motions peaks at the resonance frequency uj* 
with a characteristic speed v\, Eq. (27b). However, in 
the high-frequency dynamics we observe the co-existence 
of a small, but non-vanishing time-averaged compo- 
nent, < max r {{v(r, t)) t } <§C v\. Fig. 4(a) shows 
max r {(v(r, t)) t } as a function of frequency for different 
external voltages. The corresponding corresponding solid 
line shows exact results within the Debye-Hiickel approx- 
imation [20] . Panel (b) shows a particular example of the 
time-averaged velocity field (v(r, t)) t . 



VI. DISCUSSION 

We have analyzed the full non-equilibrium electro- 
hydrodynamics of the Debye screening layer that arise 
in an aqueous binary solution near a planar wall when 
applying a spatially modulated ac- voltage Vq cos(qy)e lu>t , 
Eq. (1). Using first order perturbation theory we have 
obtained analytic solutions for the pressure and velocity 
fields of the electrolyte and for the electric potential. Our 
analytical solution applies to the low-frequency Debye- 
Hiickel regime where the amplitude Vq of the external po- 
tential is lower than the thermal voltage Vt and the driv- 
ing frequency uj is lower than the inverse response-time 
of the electrolyte ui D = cr^/e (see Sees. IV A and IV C). 
It should be noted that our analysis does not cover the 
special case of suddenly applied step voltages, where the 
system selects its own intrinsic time scale different from 
the external time scale 1/uj [19]. 

Furthermore, we have limited ourselves to the diffusive 
regime where convection can be neglected corresponding 
to a su fficiently low driving amplitude, Vq <C V c where 
V c = y/(l + 8)rjD/e is a convective voltage (see first para- 
graph of Sec. IV B, v ~ v^y in the Debye layer). We 
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FIG. 4: Numerical simulations of flow. Panel (a) shows the maximal velocity of the harmonic flow max T ,t{v(r, t)} and the 
maximal velocity of the time-averaged flow max r {(v(r,t)}A, both versus normalized frequency. The solid lines show exact 
results within the Debye-Hiickel approximation [20] and the dashed line shows (-^ + as suggested by Eq. (38). Panel 

(b) shows an example of the time-averaged velocity field (v(r, t)) t - 



have also considered the low-frequency regime u <C uj c 
where viscosity dominates over inertia (see Sec. IV E). 

Finally, we have considered the limit with the spa- 
tial modulation being much longer than all other length 
scales in the problem, i.e., qd <C 1 and gAo < 1 (see 
Sec. IV C). In summary this means that the analytical 
studies of the effect of Eq. (1) are valid in the limits 



q <C min {d 1 , A^, 1 }, 
u <C min {oj d , uj c }, 
V < min{V T , V c }. 



(40a) 
(40b) 
(40c) 



As a main result we have supplied a proof for the validity 
of the capacitor model. The full dynamics seems however 
not to be captured by the capacitor model. Taking the 
time-average in Eqs. (27a) and (38) we get (F) 4 = 
and (vV = (in full agreement with the discussion in 
Ref. 2). In contrast, we obtain (F) ^ if we begin from 
Eqs. (18a) and (20) without expanding in oj/uj d and qXn, 
the result being finite even in the zero- frequency limit. 
Somewhat similar results were reported in another non- 
equilibrium study [5], though for a different geometry. 
Naively, this observation could suggest that ('v) t ^ 
contrary to the statement in Ref. 2. However, by also 
averaging over the y-direction we get (Fy) ty = sug- 
gesting that (Wj/) t = in agreement with the symmetry 

arguments emphasized in Ref. 2. If the finite (F) does 
not give the fluid a directional flow globally, we might 
speculate that, at high frequencies, it makes the fluid cir- 
culate in non-oscillating vortices, see Fig. 4(a), whereas 
the fluid is at rest at zero-frequency - despite (F) being 
finite. The solution to this apparent contradiction lies 
in the pressure, which will compensate the body-force at 
low frequencies. We can explicitly show that the time- 



averaged body-force can be written as 

(F(r,t)) t = Vp P (r) + 0([w/6 



(41) 



where 
PfO) = 



2-en 2 V 2 



q\ D cosh(gd) + f-y/l + {qX D ) 2 smh(qd) 



Pi 1 



x e-V'+t^) 2 ^ cos 2 (gy), (42a) 
C D \ fC cS (0)V l + (qX D ) 2 



c D 



(qX 



D 



x e -2^/i+(^oYx/\ D cos 2( (?2/ )_ (42b) 



The form of Eq. (41) suggests that 

<v(r,i)) t = + O([c/^] 2 ) 



(43) 



with p F being a pressure that compensates the low- 
frequency part of the body force, see Eq. (30). The time- 
averaged velocity field (v(r,i)) t can be calculated rig- 
orously and the complex expression (not shown) agrees 
fully with Eq. (43). At high frequencies we expect sta- 
tionary vortices, Fig. 4(b), to co-exist with the harmonic 
fluid motion illustrated in Fig. 2(c), whereas at low fre- 
quencies the circulation vanishes and we are left with 
the pure harmonic motion. Our time-dependent finite- 
element simulations in Fig. 4 support this picture and 
similar time-averaged flow in a slightly different geome- 
try has been observed both experimentally, theoretically, 
and numerically [4-6]. 



VII. CONCLUSION 

Our results provide the theoretical underpinning of 
the capacitor model widely used in the literature [2- 
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4, 6, 7, 15, 19], and form a firm starting point for fu- 
ture studies of electro-kinetic pumps and mixers driven 
by spatially modulated surface potentials. In general for 
large values of 5 we find that the Debye-Hiickel approx- 
imation works well even at elevated voltages [agreement 
within less than (1 + <5) _1 ] where it is typically expected 
to work poorly. However, our non-equilibrium approach 
has also revealed interesting short-comings in the capac- 
itor approach for high-frequency dynamics where static 
vortices appear along with the harmonic rolls also pre- 
dicted by the capacitor model. 
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